Little-Parks effect in a superconducting loop with magnetic dot 
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We have studied the nucleation of superconductivity in a mesoscopic Al loop, enclosing magnetic 
dot with the perpendicular magnetization. The superconducting phase boundary T C (B), determined 
from transport measurements, is asymmetric with respect to the polarity of an applied magnetic 
field. The maximum critical temperature has been found for a finite applied magnetic field, which is 
antiparallel to the magnetization of the dot. Theoretical phase boundary shows a good agreement 
with the experimental data. 

PACS numbers: 74.78.Na,73.23.-b,74.25.Dw 

A great deal of experimental and theoretical work has been devoted to the investigation of mesoscopic supercon- 
ductors (e.g. ref. [l|,|2| and references therein). Over the past few years hybrid superconductor/ferromagnet systems 
have attracted a lot of attention [E 0, IE IE IE IE IE EH EH E3- It has been revealed that a magnetic dot, embedded 
into a superconductor, attracts a vortex, resulting in an enhanced pinning when the magnetic moment m of magnetic 
dot and an applied magnetic field B have parallel orientation Moreover, it was recently demonstrated that 

perpendicularly magnetized magnetic dots could efficiently be used not only to enhance the flux pinning, but also to 
tune the nucleation of superconductivity |12j . 

Theoretical studies of individual mesoscopic superconductor/ferromagnet structures in the disk geometry were car- 
ried out in the framework of the Ginzburg-Landau formalism, with the emphasis on the vortex phases and transitions 
deep in the superconducting state [6J, A mesoscopic superconducting Al disk with a perpendicularly magne- 
tized magnetic dot was recently fabricated and its superconducting phase boundary T C (B) was measured Q. These 
measurements demonstrated that the superconducting critical temperature is systematically higher when the applied 
magnetic field is parallel to the magnetization of the dot. However, in this sample the superconducting phase boundary 
seems to have been affected by the proximity effects between the superconducting disk and magnetic dot. 

In this paper the nucleation of superconductivity in a mesoscopic superconducting Al loop, with a magnetic dot 
at the centre of the opening, is investigated. Given the spatial separation of the superconducting loop and magnetic 
dot, it is clear that the dot interacts with the loop exclusively through the stray field and proximity effects are 
absent. Therefore, any peculiar features in the T C (B) phase boundary can undoubtedly be attributed to the magnetic 
interaction between the superconducting condensate and the stray field of the dot. 

The sample was prepared by electron beam lithography in three steps on a Si02 substrate using positive PMMA950K 
and the corresponding co-PMMA electron beam resists. The lithography was carried out in a JEOL 5600 scanning 
electron micrograph modified for electron beam writing. The contacts pads, with the size of 500 pan, leads and 
alignment markers were patterned first and 5 nm of Cr and 30 nm of Au were thermally evaporated. Afterwards, the 
resist was removed with a lift-off procedure, using warm acetone and ultrasonic agitation. The sample was covered 
again with the e-beam resist. After the alignment, the mesoscopic loop and contacts were patterned. The designed 
width of the loop and contacts is 400 nm, whereas the inner radius of the loop is 600 nm. In order to avoid proximity 
effects, the mesoscopic contacts next to the loop are 22 /jm apart from the Cr/Au leads. 35 nm of Al was thermally 
evaporated and a lift-off was carried out. In the final step, the sample was yet again covered with the e-beam resists. 
After a quite precise alignment procedure, magnetic dot with the radius of 200 nm was patterned. A magnetic dot 
consisting of 2.5 nm Pd buffer layer and 10 bilayers of 0.4 nm Co and lnm Pd was subsequently evaporated. The 
resist was removed in a lift-off. 

Fig. ^ shows an electron micrograph of the structure. The bright area at the centre is magnetic dot. The actual 
dimensions of the loop and magnetic dot were obtained from atomic force microscopy. The inner radius of the loop 
equals Ri = 0.55 /im, outer radius is equal to R = 1.05 ^m, whereas the Co/Pd magnetic dot has the radius of 
rd = 0.27 /im. The mesoscopic contacts next to the loop, used for measuring voltage, are closer to the loop than 
initially designed. This effectively increases the outer radius of the loop, since the supercurrent can spread into the 
contacts, thereby affecting both the periodicity and background of the phase boundary Q]. 

The superconducting phase boundary T C (B) was obtained by four-point transport measurements in a cryogenic 
setup at temperatures down to 1.11 K, applying the magnetic field perpendicularly to the sample surface. A transport 
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FIG. 1: A scanning electron micrograph of the structure. 



current with the rms of 150 nA and frequency 27.7 Hz was used. The response was measured by a PAR124A lock- 
in amplifier. The phase boundary was determined resistively by sweeping the magnetic field at a very low rate, 
whilst keeping the temperature constant. T C (B) curve was afterwards extracted from the magnetoresistance using 
the resistive criterion R n /2, where R n is the resistance in the normal state. The temperature and field steps were 
500 fiK and 5 /iT, respectively, with the temperature stability 100 fiK. In order to minimize the influence of possible 
temperature gradients, the temperature was being monitored by two independent temperature sensors placed at two 
different positions on the sample holder. 

The resistance of the sample at room temperature is 6.6 Q, the resistance in the normal state at low temperatures 
is R n = 2.5 17, whereas the maximum critical temperature T cm is 1.3422 K. 

The multilayer Co/Pd magnetic dot provides a perpendicular magnetization |l2l Il4|. The hysteresis loop in the 
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FIG. 2: The hysteresis loop of a co-evaporated Co/Pd plane film obtained by the magneto-optical Kerr measurements at room 
temperature. The shaded area indicates the range of the applied magnetic field. 

perpendicular magnetic field of the co-evaporated plane film with the same composition, obtained by the magneto- 
optical Kerr measurements at room temperature, is presented in Fig. [21 The shaded area indicates the range in 
which the applied magnetic field was varied in the experiment. The remanence is complete, whereas the coercive 
field is approximately 70 mT. When saturated, Co/Pd structures patterned at the sub-micrometre scale are in the 
single-domain state [l2j . 

Fig. shows the calculated stray field produced by a Co/Pd magnetic dot, averaged over the thickness of the 
superconducting loop. The stray field is normalized to the magnetization of the magnetic dot. The darker shaded area 
shows the position of the magnetic dot, whereas the lighter shaded areas indicate the position of the superconducting 
loop. The arrow shows the direction of the magnetization m of the dot. The stray field was calculated assuming 
a single-domain state of the magnetic dot. Prior to the measurements the dot was magnetized perpendicularly in 
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FIG. 3: The calculated stray field of the magnetic dot, averaged over the thickness of the superconducting loop. The stray field 
is normalized to the magnetization of the magnetic dot and was calculated assuming a single domain structure. The darker 
shaded area indicates the position of the dot, whereas the lighter shaded area shows the position of the superconducting loop. 
The arrow shows the direction of the magnetization. 



300 mT. As the applied magnetic field in our experiment never exceeded 30 mT (Fig. |2J, it has been assumed that 
the magnetization of the magnetic dot remains unaltered during the measurements [l2( • 

The fluxes of the stray field (see Fig. |3J) through the opening of the loop (r < 0.55 /xm) and loop itself (0.55 < 
r < 1.05 [/jm]) have different orientations and, in qualitative terms, there should be a competition between the 
effects of these two fluxes on the superconducting condensate. Assuming that the loop is sufficiently big compared 
to the coherence length at zero temperature £(0) ( Ri,R >> £(0), where Ri and R are the inner and outer radius 
of the loop, respectively) the loop is similar to a stripe at temperatures much lower than the maximum critical 
temperature the system can attain T cm (we will be referring to these temperatures as low temperatures). Clearly, 
the superconducting state of the stripe depends only upon the flux value within the stripe. In other words, at low 
temperatures the behaviour of the loop is governed by the flux through the loop and not by the flux in the opening 
of the loop. Therefore, at low temperatures the shift of the T C (B) phase boundary along the £?-axis is expected to 
be in the direction of the dot magnetization. At temperatures close to the maximum critical temperature (we will 
be referring to these temperatures as high temperatures), the direction of the shift of the T C {B) phase boundary is a 
result of an interplay between the opposite fluxes. Note that the superconducting phase boundary does not depend 
on the local profile of the stray field in the opening of the loop, but on the total flux through the opening generated 
by magnetic dot. 

To describe quantitatively the T C {E) phase boundary the Ginzburg-Landau theory has been used. The free energy 
of the superconducting structure can be expressed as 
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where tp is the modulus of the order parameter, A is the vector potential, an integer number L stands the winding 
number (vorticity) of the order parameter, $o is the superconducting flux quantum, whereas a and p represent the 
Ginzburg-Landau parameters. The integration is performed over the area of the loop. In Eq. cylindrical symmetry 
of the order parameter near the phase transition to the normal state has been taken into account. As the thickness 
of sample is much smaller than £(T) within the temperature range of interest, there is no modulation of the order 
parameter in the direction of the applied field. For this reason, the stray field can be averaged out over the thickness of 
the sample and the problem is reduced to the ID case. Near the superconducting phase boundary the total magnetic 
field is equal to the sum of the stray field and an applied magnetic field. Therefore, the vector potential may be 
expressed as A(r) = A s (r) + Br/2, where A s (r) is the vector potential of the stray field and B is an applied magnetic 
field. We represent ip(r) as a Fourier expansion within the range Ri < r < R Q and consider the coefficients of the 
expansion the variational parameters. The values of these parameters can be calculated by a minimization of the free 
energy, yielding the solution of the full Ginzburg-Landau system of equations. The transition to the normal state 
corresponds to ip = 0. We have found that for our sample and experimental conditions, only the first term in the 
Fourier expansion is important, that is, tp is constant within the sample, whereas all other terms are negligibly small 
and practically do not affect the position of the phase boundary. 
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FIG. 4: The measured (open symbols) and calculated (solid line) superconducting phase boundary given as 1 — T C (B)/T cm 
versus the normalized flux &/&o ■ Here T cm is the maximum critical temperature and $o is the superconducting flux quantum. 



Fig. 01 presents the superconducting phase boundary, given as 1 — T C (B)/T cm versus the normalized flux $/$o- 
T cm is the maximum critical temperature and <i>o is the superconducting flux quantum. Applied fields parallel to the 
magnetization of the dot are taken as positive and vice versa. The open symbols are experimental data, whereas the 
solid line corresponds to the theoretical curve for Ri = 0.55 /im, and rd = 0.27 fj/m. The coherence length at zero 
temperature £(T = 0) has been used as a fitting parameter as it depends upon the preparation procedure and can 
show moderate variations even amongst samples prepared simultaneously. Due to the influence of the mesoscopic 
contacts, the outer radius has to be used as a fitting parameter, since the supercurrent is not confined only by the 
ring, but can well flow into the mesoscopic contacts, thus effectively increasing the outer radius. The best fit has been 
obtained for the coherence length at zero temperature £(0) = 96 nm and the outer radius R = 1130 nm. The latter 
is quite close to the actual radius R D — 1050 nm. 

The theoretical T C {B) curve is in a good agreement with the experimental data. The asymmetry, shift and periodicity 
of the experimental data have been reproduced well. The asymmetry and the shift can be explained by considering 
Eq. (I). From this equation it follows that the maximum critical temperature, for a fixed L, is attained when the 
applied field is equal to 
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The corresponding critical temperature is 
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where T coo stands for the bulk critical temperature of Al. For the sake of simplicity, we assume that the stray field 
does not change considerably over the width of the loop. Denoting the average field in the opening of the loop and 
within the loop as B op and B t , respectively Eqs. |J2J and J2J reduce to 
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The maximum critical temperature T cm can be found by the minimization of the right-hand side of the Eq. JSJ with 
respect to the integer L. Eq. Q gives the shift of the T C (B) curve along the B-axis close to the T cm , caused by the 
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presence of an additional magnetic field. This shift is a sum of two contributions. The first term is just a compensation 
of the stray field B[ within the loop. The second term describes the effect of the inhomogeneity of the stay field. This 
contribution behaves rather nontrivially due to the discreteness of L. For instance, when < ttRJ{B op — B{) < 0.5$o, 
it follows from Eq. (5) that L — 0, and the second contribution to the shift is negative. If the size of the dot and/or 
the intensity of the stray field are increased 0.5$o < irB%(B op ~ Bi) < <I>o, L switches to 1, and the second contribution 
becomes positive. Physically, this means that it is energetically more favorable for the system to change the vorticity 
by f in order to decrease the modulus of the current in the loop. As a result, the current and the total shift, defined 
by Eq. (4), can change the sign. Thus, we may conclude, that the contribution to the shift, which is caused by the 
stray field inhomogeneity, is a periodic function of B op — Bi and vanishes in the case of a homogeneous field, B op = Bi . 
The total shift of the T C (B) curve along the B-axis close to the T cm can be of both signs for the same orientation of 
magnetization of the dot. 

Our calculations have shown that at low temperatures the shift of the phase boundary is controlled only by the 
average stray field within the loop and is equal to -Bi, which is in accordance with the preliminary qualitative 
analysis. If the additional magnetic field were homogeneous, the shifts of the T C (B) phase boundary for low and 
high temperatures would have the same direction. The inhomogeneity of the stay field gives rise to two different and 
mutually independent shifts of the phase boundary in the low and high temperature ranges, which, in turn, twist the 
whole superconducting phase boundary, as presented in Fig. 0] 

We would like to note that the discrepancy between the theoretical and the experimental data for low temperatures 
can be attributed to the effect of contacts, since they change the geometry of the sample, which cannot be entirely 
accounted for only by using the outer radius of the loop as a fitting parameter. 

To summarize, we have fabricated a mesoscopic superconducting loop with a magnetic dot at the centre of the 
opening and investigated the onset of superconductivity in this structure by measuring the superconducting T C (B) 
phase boundary. The phase boundary has been found to be asymmetric with respect to the polarity of an applied 
magnetic field. The structure exhibits the maximum critical temperature for a finite value of the applied magnetic 
field, which is antiparallel to the magnetization of the magnetic dot. The theoretical superconducting phase boundary, 
obtained in the framework of the Ginzburg-Landau theory, is in a good agreement with the experimental data. It has 
been demonstrated that the inhomogeneous stray field shifts the T C (B) phase boundary along the B-axis independently 
in the low and high temperature ranges. In the limit of a thin loop with a large radius, these shifts are determined 
by the average value of the stray field in the opening of the loop and within the loop. 
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